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5^ The excited states of 12 Be, 14 Be and 15 B were studied by an antisymmetrized 

Q-i. molecular dynamics method. The theoretical results reproduced the energy 

levels of recently measured excited states of 12 Be, and also predicted rotational 
bands with innovative clustering structures in 12 Be, 14 Be and 15 B. Clustering 
states with new exotic clusters ( 6 He, 8 He and 9 Li) were theoretically sug- 
gested. One new aspect in very neutron-rich nuclei is a 6-nucleon correlation 
among 4 neutrons and 2 protons, which plays an important role in the for- 
\Q , mation of 6 He clusters during clustering: 8 He+ 6 He of 14 Be and 9 Li+ 6 He of 

15 B. 

^ Owing to progress of experimental techniques, information concerning the excited states 

of light unstable nuclei has increased rapidly. Clustering structures should be one of the 
essential features in unstable nuclei as well as in stable nuclei. However, the characteristics 
of clustering in unstable nuclei are still mysterious. Do clustering structures often appear in 
the excited states of general light unstable nuclei ? Searching for exotic clustering structures 
and other new-type clusters than the well-known a cluster is the focus in studies of the 
excited states of light unstable nuclei. The clustering structures of unstable nuclei are one 
of the attractive and important subjects in experimental and theoretical research. 

Theoretical studies on 10 Be [|l],0||,f|,|5| have suggested molecule-like structures in the 
K = band and in the excited K = O2" band. Other candidates of clustering structures 
are the excited states of 12 Be observed recently in the inelastic scattering reactions || and 
in break-up reactions into 6 He+ 6 He and 8 He+ 4 He channels [[F]J|. However, fully microscopic 
calculations have not yet been done for these states of 12 Be. Another interesting feature of 
12 Be is a vanishing of the neutron magic number N = 8 PJlQ]. If the ground state is an 
intruder state, where are the normal states with closed neutron-shell configurations? It is 
important to solve the relation between the intruder ground state and the possible clustering 
states within one microscopic framework in order to systematically investigate the structures 
of 12 Be. Our first aim is to study the ground and excited states of 12 Be with a microscopic 
method. Furthermore, it is a challenge to theoretically search light neutron-rich nuclei ( 14 Be 
and 15 B) for new clustering states with exotic clusters 

We have applied a microscopic method of antisymmetrized molecular dynamics (AMD). 
The AMD methods have already proved to be a powerful approach for structure studies of 
general light nuclei 11,11,12,13] because of the flexibility of the wave functions, which can 



express various kinds of structures from spherical shell-model-like structures to developed 
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clustering structures. It is possible to systematically investigate the ground and excited 
states of light nuclei with a new version: a variation after spin-parity projections (VAP) 
of AMD, which was proposed by the author fll^fl. As for theoretical predictions with this 



method, the predicted 5~ and 6 + states of 10 Be with developed clustering structures || have 
actually been discovered in recent experiments ||. 

In the present work, we have studied the structures of the ground and excited states 
of 12 Be and compared the theoretical results with the experimental data. We theoretically 
searched 14 Be and 15 B for unknown excited states with new clustering structures. In the 
results for 14 Be and 15 B, the prediction of new rotational bands of exotic clustering structures 
are presented. We discuss the mechanism of clustering development and the origin of the 
formation of 6 He clusters. 

An AMD wave function is a Slater determinant of Gaussian wave packets: 
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where the centers of Gaussians Zj's are complex variational parameters. Xi is an intrinsic spin 
function parameterized by which is also a variational parameter specifying the direction 
of the intrinsic spin of the ith nucleon. is an isospin function which was fixed to be up 
(proton) or down (neutron) in the present calculations. We varied the parameters Zj and 
= 1 ~ A) to minimize the energy expectation value for the spin-parity eigenstate (VAP 
calculations), 
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where the operator of total-angular-momentum projection (P i 
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The integration for the Euler angle fl was calculated numerically, 
cooling method 



) is fdQD^ Kf (Q)R(n). 
We adopted a frictional 
we represent the intrinsic 



TTJ to obtain the minimum energy states. Here, 
AMD wave function as < ^>q ± , that was obtained by a VAP calculation for the lowest state with 
a given spin parity of J^. Higher excited states were constructed by varying the energy of 
the orthogonal components to the lower states by superposing wave functions. More details 
of the AMD method for the excited states with the variation after spin-parity projection 
are described in Refs |5|,[n|]. By performing VAP calculations for the lowest and the higher 
excited states with various sets of total spin and parity {J±}, we obtained many AMD 
wave functions {$i, • ■ ■ , $ m }, each of which approximately represents the intrinsic state of 
the corresponding state. The final results were attained by diagonalizing a Hamiltonian 



matrix, (PMK'®i\ H \ P MK»®j) &3 = 1 

We applied the AMD method for the excited states of 12 Be, 14 Be and 15 B. The adopted 
interactions in the present work were the central force of the modified Volkov No.l with case 
3 [13], the spin-orbit force of G3RS jnj and the Coulomb force. The Majorana parameter 
used here was m = 0.65, and the strength of G3RS force was chosen to be U\ = —u 2 = 3700 
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MeV. We choose an optimum width parameter (y) for the Gaussians of the single-particle 
wave functions of each nucleus. 

The theoretical results of the binding energies of 12 Be, 14 Be and 15 B are 61.9 MeV, 
59.7 MeV and 73.1 MeV, which underestimate the experimental values: 68.65 MeV, 69.77 
MeV, 88.19 MeV, respectively. We find that the smaller values of the binding energies in 
the present results than the experimental data can be improved by changing the Majorana 
parameter of the interaction small to be m = 0.61. With this parameter, the binding energy 
of 14 Be is 67.7 MeV. In spite of the smaller values of binding energies, we adopt the parameter 
m = 0.65 in the present work, because this parameter gives a good reproduction of parity 
inversion of 11 Be, which should be important to well describe the ground-state properties 
of neighboring nucleus: 12 Be. We have checked that the change of m parameter has no 
significant effect on the excitation energy -E(O^) of 14 Be, at least. Namely, the calculated 
values of E{0f ) i n the cases m = 0.61 and m = 0.65 are 4.2 MeV and 4.4 MeV, respectively. 

In the calculated results for 12 Be with AMD, many excited states appear in the low-energy 
region. The energy levels of 12 Be are presented in Fig. |]. The theoretical levels of the 4f, 62" 
and Sf states correspond well to the recently observed excited states 0. By analyzing the 
intrinsic AMD wave functions, we can classify the excited states into rotational bands, such 
as K n = Of, Of, O3", lj". It is surprising that the newly observed levels ]7) at an energy region 
above 10 MeV belong to the third rotational band, K = Of , which is a 2hu excited state 
with a well-developed clustering structure. On the other hand, there are many low-lying 
positive-parity states which belong to the lower rotational bands, K 71 = Of and K n = Of- 
Even though 12 Be has a neutron magic number of 8, the intrinsic state of the ground + 
state is not an ordinary state with a closed neutron p-shell, but a prolately deformed state 
with a developed clustering structure. The deformed ground state is dominated by another 
2p — 2h configuration than that in the 0% state. As a result, the ground K n = + band starts 
from the ground + state and reaches the band terminal at the 8f state. The rotational 
bands, K ± = Of and K ± = Of, share the same band terminal state because the terminal 
J ± = 8f state is the highest spin state in the 2hu configurations. On the other hand, the 
main components of the 0^ and 2 + states are the Ohu configuration with a closed neutron 
p-shell, and construct the second K n = Of band. 

For experimental evidence, the strength of the f3 decay from 12 Be(0^) to 12 B(1 + ) is helpful 
to estimate the breaking of the neutron p-shell [[|. As expected, the theoretical value of 
B(GT) = 0.8 is sufficiently small as the experimental data, B(GT) = 0.59, because the 
component of the 2tkv configurations in the parent 12 Be(0i~) makes the transition matrix of 
the Gamov- Teller operator to be small. Such a result is consistent with pioneer work |9,16 



where their treatments were not fully microscopic for the freedom of all nucleons. The 
calculated value of B(E1) from l^f to 0^ is 0.02 e 2 fm 2 , which agrees reasonably with the 
experimental large value, B(E1) = 0.05 e 2 fm 2 [fLOR . As for the E2 transition strength, the 
predicted values in the present results are B{E2; 2f — > 0^) = 14 e 2 fm 4 and B{E2; 2f — * 0^) 
= 8 e 2 fm 4 , which indicate a larger value of the E2 transition strength in the deformed 
K = Of band than that in the K = Of band. 

It should be pointed out that this is the first theoretical work which systematically 
reproduces the energy levels of all spin-assigned states discovered recently, except for the 1~ 
state. Although the excitation energy of the 1~ state is overestimated by theory, we regard 
the theoretical lj~ state as being the measured 1~ state at 2.68 MeV because of the large 
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FIG. 1. Excitation energies of the levels in 12 Be. The theoretical results were calculated by 
VAP calculations based on AMD. The experimental data are from the Table of Isotopes and Refs 



B(E1) value. A theoretical discovery in the present calculations is the K n = Of rotational 
band with the developed 6 He+ 6 He clustering structure, which constructs the Of, 23 , Af, 62" 
and 8f states. The developed clustering in the O3" state gradually weakens with increasing 
total-angular momenta (J), and finally it changes to a spin-aligned state at J" =l= = 8^. Many 
excited states in the rotational bands K 71 = Of, Of and 1~ are predicted to exist in the 
low-energy region with an excitation energy below 10 MeV. 

By analyzing the single-particle wave functions in the intrinsic state, the deformation 
mechanism of the ground state can be understood by the idea of molecular orbits surrounding 
2a cores. Although the molecular orbits in Be isotopes have been studied theoretically 
|T|,PIP],|r7f| , one of the important points of the present results is that the development of a 
2a+4n structure in the ground state of 12 Be can be confirmed without assuming the existence 
of clusters. On the other hand, the third + state seems to be a 6 He+ 6 He clustering structure 
instead of the a and surrounding neutrons in molecular orbits, because the distance between 
two clusters is too large. In this case, an a cluster goes outward far enough to form a 6 He 
cluster with 2 valence neutrons. 

We next present the results for 14 Be. Figure shows the energy levels of positive-parity 
states of 14 Be obtained by VAP calculations. Although few excited states are experimentally 
known, many excited states are predicted based on the theoretical results. There are a few 
candidates of levels above the threshold energy (9.1 MeV excitation) for separation into 
8 He+ 6 He in break-up experiments 



The excited states in the region J < 6 are classified into 3 rotational bands: K = 0f, 
K = 0~2 and K = 2 + . In most of the states in 14 Be, 2-a cores are formed in the calculated 
results, although no clusters are assumed in the model. Because of the 2-a cores, the Of 
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and Of states have prolate deformations with deformation parameters of (3 = 0.49 and 
(3 = 0.64, respectively. The latter one with a large deformation constructs the K n = Of 
rotational band, which reaches the state. By analyzing the single-particle wave functions 
of the intrinsic states, it is found that the states in the K = Of and K = 2 + bands are 
dominated by the normal OTiuj configurations, while the largely deformed states in K w = Of 
mainly contain 4p-2h states with 2hu configurations. The unique result of 14 Be is innovative 
clustering structures of the excited states in this second K 77 = Of band. The intrinsic states 
of K n = Of band have spatially developed 8 He+ 6 He clustering state, where both clusters 
are, themselves, very neutron-rich nuclei. In Figure |^(a) for the contour surface for the 
matter density, p > 0.16 nucleons/fm 3 , of the intrinsic state of the 0^ state, the 8 He( 6 He) 
cluster can be clearly recognized in the right (left)-handed side. 

In order to understand the origin of 8 He+ 6 He clustering formation, it is valuable to in- 
spect the single-particle behavior of the valence neutrons. As mentioned above, the 8 He+ 6 He 
clustering state in the K n = Of band comes from the 4p-2h configurations. According to 
analyses of the single-particle energies and wave functions, the highest 4 neutron orbits cor- 
respond to two kinds of spatial orbits, which are associated with the s<i-mixing orbits in the 
deformed system. Roughly speaking, each of two kinds of spatial orbits is occupied by a 
spin-up neutron and a spin-down neutron. Figures |3](b) and |3|(c) present the higher s<i-like 
orbit and the lower s<i-like orbit, respectively. They contain more than 80% components of 
the positive parity states. If we call the longitudinal direction of the prolate deformation 
as the z-axis, the highest nucleon orbit seen in Fig. 0(b) originates from the s(i-orbit with 
a form of yzexp[— vr 2 \. On the other hand, the spatial orbit for the lower one is similar to 
the orbit of z 2 exp[— z/r 2 ], which deforms along the z-axis, as shown in Fig.|3](c). The original 
sci-orbits are modified because of the deformation in 14 Be. The deformed neutron orbits 
are stabilized by 2 stretched pairs of protons, because the orbits gain their kinetic energies 
with respect to the nodes along the z-axis. From the viewpoint of single-particle orbits, the 
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FIG. 3. Intrinsic structure of 0^ of 14 Be. Figure (a) shows a surface cut of the matter density 
p > 0.16 nucleons/fm 3 . The density was calculated for the intrinsic AMD wave function before 
a spin-parity projection. The density distributions for the single-particle orbits of the valence 
neutrons are presented in figures (b) and (c). 

6 He cluster in the 6 He+ 8 He clustering system consists of 4 neutrons in the sc?-orbits and 2 
protons in the p-orbits. In other words, the formation of 6 He cluster originates from the cor- 
relation of 4 neutrons in the so?-orbits and 2 protons in the p-orbits. Once a 6 He cluster takes 
shape, the cluster spatially develops out of the core-cluster 8 He keeping the stability of the 
modified sd-mixing orbits in the deformed system. Unique points of this idea are as follows. 
Firstly, correlations of more than 4 nucleons have not been found. The 6 He cluster is related 
to 6-nucleon correlations (4 neutrons with 2 protons). Second, the nucleons, which compose 
a cluster out of the core nucleus, rise over 2 shells (sd-shell and p-shell in this case). This is 
unusual in light stable nuclei. The third point is that the spatially clustering development 
is associated with the stability of the single-particle orbits of the system. We suggest one of 
the novel feature that 4 neutrons with 2 protons correlate to form a developed 6 He cluster 
in neutron-rich nuclei. 

In the results concerning the excited states of 15 B, the situation is similar to that of 14 Be. 
We find an exotic clustering structure with 9 Li+ 6 He in the second K n = 3/2^ band, which 
starts from the J = 3/2^ state at about 7 MeV higher than the ground state. The intrinsic 
structure of the K n = 3/2^ band consists of two exotic clusters: 9 Li and 6 He clusters, which 
develop as remarkably as those in the K n = band of 14 Be. These excited states with 
9 Li+ 6 He clustering are based on 2hu configurations. The spatial parts of 4 valence neutrons 
in sd-like orbits are very similar to those of 14 Be presented in Figs. 0(b) and 0(c). It is an 
interesting problem whether or not the 6 He-cluster may appear in the excited states of other 
neutron-rich nuclei than 14 Be and 15 B. 
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Although particle decay widths are important for the stability of the excited states, the 
widths were not considered in the present calculations, which were performed within a bound 
state approximation. The widths for the particle decays should be carefully discussed using 
other frameworks, such as a method with a reduced width amplitude, or a complex scaling 
method beyond the present AMD framework. We just mention the partial widths of the 
states into He channels in the K n = 0^ of 14 Be. When the excitation energies are above the 
threshold energy, 12.2 MeV, the excited states in K n = of 14 Be may decay into 8 He+ 6 He 
channel. By assuming 8 He and 6 He clusters in SU(3) limits, we estimated the partial decay 
width of Qt and 8^ into 8 He(0 + )+ 6 He(0 + ) channel to be 6 keV and 50 keV, respectively, 
by calculating the reduced width amplitude. For the total width, it is required to carefully 
study the excitation energies and other decay channels, such as neutron decays and excited 
He cluster decays. 

In summary, we studied the structure of the excited states of 12 Be, 14 Be and 15 B based 
on the framework of the AMD method. This is the first microscopic calculation to system- 
atically reproduce the energy levels of all the spin-assigned states in 12 Be, except for the 1~ 
state. One of the discoveries of the present results is the third K n = + band with 6 He+ 6 He 
clusters, which corresponds well to recently observed states in the break-up reactions. In 
14 Be and 15 Be, we predicted low-lying excited states with innovative clustering structures. 
Namely, new types of clusters, 8 He+ 6 He and 9 Li+ 6 He, were found in the excited states of 
14 Be and 15 Be, respectively. In these well-developed clustering structures with 2%uj excited 
configurations, the mechanism of the 6 He cluster development was discussed in relation to 
the single-particle picture. It was found that the 6-nucleon correlation among 4 neutrons 
and 2 protons plays an important role in the formation of 6 He cluster. The results suggest 
an important feature, that an exotic clustering structure may very often exist in the excited 
states of neutron-rich nuclei. 
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